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, Diffractive heavy vector meson photoproduction accompanied by proton dissociation is studied 

for arbitrary momentum transfer. The process is described by the non-forward BFKL equation, 
for which a complete analytical solution is found, giving the scattering amplitude. The impact of 
non-leading corrections to the BFKL equation is also analysed. Results are compared to the HERA 
qj ' data on J / ip production. 

1 Introduction 

Quantum Chromodynamics offers unique opportunities to study the richness of dynamical phenomena 
of nonlinear quantum field theory. One of the most interesting problems is related to the colour flow 
in high energy scattering. In particular, diffractive processes correspond to an exchange of a colour 
singlet system of quarks and gluons between scattering objects. Such diffractive phenomena possess a 
very clean experimental signature, namely a large rapidity interval devoid of particles (i.e. a rapidity 
gap). 

The perturbative QCD description of the hard colour singlet exchange across a large rapidity 
interval y relies on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation jL], Q]. In this framework, 
the leading powers of rapidity in the perturbative expansion are resummed, giving the amplitude for 
hard pomeron exchange. The pomeron is viewed as a composite system of two reggeized gluons in the 
colour singlet state. The status of the BFKL approach to QCD amplitudes is still under discussion, 
and both theoretical improvements and experimental tests are necessary. 




Figure 1: Feynman diagram illustrating vector meson photoproduction at high momentum transfer. 



Diffractive photoproduction of a heavy vector meson, separated from the proton remnant by a 
large rapidity gap has been proposed |3], ||] as an ideal probe of the BFKL pomeron, see fig. ||. Indeed, 
this process permits detailed studies of both the momentum transfer and the rapidity dependence of 
the scattering amplitude. The vector meson mass and the momentum transfer t provide the hard scale 
required for the perturbative treatment of QCD processes, and the sensitivity to the infra-red region 
is small, contrary to the case of inclusive hard diffraction. 

There are some very recent measurements of this process from HERA || that allow the theoretical 
models to be tested. The available calculations (3|, ||, |6| of the cross-sections for heavy vector meson 
production are based on the Mueller- Tang approximation Jjj] to the solution of the leading logarithmic 
BFKL equation. In this approximation, parts of the amplitude which vanish in the limit y — > oo are 
neglected [||, [§. 

Recently, it has been shown 0] that fitting these results to the data, one obtains a good quantitative 
agreement with the differential cross-section. Still, the important non-leading corrections to the BFKL 
kernel ||] are not accounted for. Furthermore, the Mueller- Tang approximation is only good for very 
large rapidities, and may need improvement in order to understand the experimental data. Indeed, 
it has been found that to describe the events with gaps between jets subleading corrections to the 
Mueller- Tang picture are important [Q, 10|. 

Thus, the main goal of this paper is to investigate diffractive heavyQ vector meson photoproduction 
beyond the leading logarithmic BFKL equation and beyond the Mueller- Tang approximation. The 
obtained results are compared to previous ones ||, |[] and to the experimental data from HERA |J. 
In Sec. |2] we define the framework, in Sec. [| the BFKL equation is presented, and in Sec. || an exact 
solution of the equation is derived. Properties of the exact and numerical solutions are studied Sec. 
||, comparison with data is performed in Sec. |6|, and in Sec. |?j conclusions are given. 
1 Light vector meson production will be studied in a forthcoming paper Jill]. 
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2 Hard colour singlet exchange 



The diffractive process 7p — > at large momentum transfer t (see fig. [l]) takes place by exchange of 
the BFKL pomeron. It has been demonstrated, that at large momentum transfer, the hard pomeron 
couples predominantly to individual partons in the proton [12]. Thus, the cross-section may be 
factorized into a product of the parton level cross-section and the parton distribution functions, 

where N c = 3, G(xj,t) and qj(xj,t) are the gluon and quark distribution functions respectively, and 
W 2 is the 7p centre-of-mass energy squared. The struck parton in the proton, initiating a jet in the 
proton hemisphere, carries the fraction Xj of the longitudinal momentum of the incoming proton. The 
partonic cross-section, characterized by the invariant collision energy squared s = XjW 2 is expressed 
in terms of the amplitude A(s, t), 

The amplitude is dominated by its imaginary part, which we shall parametrize, as in ||, |], by a 
dimensionless quantity T 

Imi(s,t) = ^(z,r) (3) 



where z and r are defined by 



(4) 

(5) 



where My is the mass of the vector meson, Q 7 is the photon virtuality^ and A 2 is a characteristic mass 
scale related to My and \t\. Following the results of Q we assume A 2 = My+Q 2 . For completeness, we 
give the cross-section expressed in terms of !F(z, r), where the real part of the amplitude is neglected, 

da(jq -> Vq) 16vr 2 

It = 81^ ( *' r)l ' (6) 

This representation is rather convenient for the calculations performed in Sec. [§. 



3 The BFKL equation 

The imaginary part ImA(s, t) of the amplitude for the process ■jp — ► y + gap + X+jet corresponds to 
the diagram in fig. illustrating QCD pomeron exchange, and can be written in the following form: 

2 . _ fd 2 k ^ v (k 2 ,q 2 )^ qq (x,k,q) 



2m A(s,t = -q)=l — [{k + q/2? + sm _ q/2? + SQ] - W 



2 In this paper we only consider Q 1 = 0. 
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In this equation, x is the longitudinal momentum fraction of the incoming proton taken by the hard 
pomeron, x = A 2 /s, q/2 zk k denote the transverse momenta of the exchanged gluons, and q is 
the transverse part of the momentum transfer. In the propagators corresponding to the exchanged 
gluons we include the parameter sq which can be viewed as the effective representation of the inverse 
of the colour confinement radius squared. Investigation of such non-perturbative effects in the gluon 
propagator at low virtualities was performed, for instance, in [13], where it was found that 0.1 GeV 2 < 



sq < 0.5 GeV . Recent lattice studies [14| indicate a value of sq between 0.25 GeV and 0.65 GeV . 



Having some freedom here, we choose sq = 0.5 GeV 2 . The sensitivity of the cross-section to the 
magnitude of so can serve as an estimate of the sensitivity of the results to the contribution coming from 
the infra-red region. It should be noted that formula (0) gives a finite result in the limit sq = GeV 2 . 
An interesting discussion of non-perturbative effects in the gluon propagator in the context of vector 
meson production can be found in |]15||. 



The couplings of the external particle pair to the colour singlet gluonic ladder are described, in 
the high energy limit, by impact factors ^yik 2 ,q 2 ) and &® q (k 2 , q 2 ) for the 7 — > V transition and 
the quark elastic scattering, respectively. The impact factors are obtained in the perturbative QCD 
framework and we approximate them by the leading terms in the perturbative expansion [p"G|]: 

$ o Ca * ( 1 1 



l y 2 \q 2 q^ + k 2 , 

*° qq = (8) 

In the former formula, factorization of the scattering process and the meson formation is assumed, 
and the non-relativistic approximation of the meson wave function is used. In this approximation 
the quarks in the meson have collinear four-momenta and My = 2M q where M q is the mass of the 
constituent quark. To leading order accuracy, the constant C may be related to the vector meson 
leptonic decay width 

C 2 = (9) 

We have also defined 

f = qj + q 2 /4, (10) 
q 2 = (Q 2 +M 2 )/4. (11) 

This standard approximation is based on the fact, that the typical quark three- momenta in charmo- 
nium states are much smaller than the charm quark mass. Still, these momenta are non-zero, leading 
to a smearing of the non-relativistic form- factor. The resulting corrections were extensively discussed 
e.g. in [[l?]] for the forward scattering case. Moreover, it is expected at larger momentum transfer 
\t\ > My, that the relativistic effects become enhanced, as the mass is no longer the largest scale in 
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the process. For instance, in the non-relativistic approximation one neglects helicity flip processes, 
whose relative contribution may grow with increasing momentum transfer, as was demonstrated for a 



single helicity flip amplitude for light vector mesons [18|. Therefore, the applied approximation may 
be inaccurate at very large momentum transfers. 

The function & qq (x, k, q) satisfies the BFKL equation, which in the leading ln(l/x) approximation 
has the following form: 

d 2 k' 



2 (k'-k) 2 + S 



kf + s k'i + s 

.11 



(k? + s )(k% + so 



12 



+ 



k\ 

72 



k'{ + (kf - k) 2 + 2s k'i + (k' - k) 2 + 2s, 



- k) 2 + s 

<f>qq(x',k',q)- 

$ qq (x',k,q)\ 



where 



=1,2 



?±fc. 



«±h> 

2 



(12) 



(13) 



denote the transverse momenta of the gluons. At leading logarithmic accuracy, a fixed value of the 
QCD coupling a s should be used in equations (||) and (p^). 

It is known that the BFKL equation can acquire significant non-leading contributions |J . Although 
the structure of those corrections is fairly complicated, their dominant part is rather simple, and follows 
from restricting the integration region in the real emission term in equation (^) [jl^, pQ ]. For q = 
the relevant limitation is [19, ^] 

(14) 



k' 2 < k 2X -. 



This follows from the requirement that the virtuality of the gluons exchanged along the chain is 
dominated by the transverse momentum squared. The constraint (|l4|) can be shown to exhaust about 
70% of the next-to-leading corrections to the QCD pomeron intercept H, Q9|. Generalization of the 
constraint (|14|) to the case of a non-forward configuration with q 2 > is assumed to take the following 
form ||, 0: 



k < (k 2 + q 2 /4)- 



x 



(15) 



The latter formula gives at k 2 < g 2 /4 a less restrictive relation for k' in respect to k, than constraint 
( |l4|) in the forward case. Consequently, the generalized constraint (^) may realize, at larger momen- 
tum transfer, a smaller fraction of the non-leading corrections than the quoted 70% at q 2 = 0. This 
would mean that our estimate of the cross-section could be too high at larger q 2 . 

Another important part of the non-leading corrections to the BFKL equation is related to running 
of the coupling constant within the ladder. To be consistent, the running coupling will also be used 
in the impact factors (g). Besides the BFKL equation (|i~2| ) in the leading logarithmic approximation 
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we shall therefore also consider the equation embodying the constraint ( |15| ) and a running coupling 
in order to estimate the effects of the non-leading contributions. 



The corresponding equation which contains constraint ( |15| ) in the real emission term reads: 

3a s (n 2 ) f 1 dx' f d 2 k' 



<S> qq (x,k,q) = <S> q ] (k 2 ,q 2 ) + 



1 d^_ 

x' 



k\ 



+ 



kf + s k'i + s 



2 jk' - k) 2 + sp 

^ "J2 , „ \/i„/2 



(k[ 2 + s )(k' 2 + s 



2vr 2 J x x' J (k'-k) 2 + s 

% q (x', k', q)G {{k 2 + q 2 /A)x'/x - k' 2 )) 



k\ 



+ 



k\ 



(16) 



_kf + (k' - fc) 2 + 2s k'i + (k' - k) 2 + 2s . 
with the scale of the coupling set to fi 2 = k 2 + q 2 /4 + so- The scales of the coupling constants in the 
impact factors should be related to the virtualities entering the vertices. A natural choice is then 



/i 2 = k 2 + M 2 



m 



/j,?, = k 2 + s 



111 



qqi 



(17) 



with k 2 being the virtuality of the gluon entering the vertex. We will also consider another choice 



Mi 



/xf/4 and // 2 = //|/4, 



(18) 



necessary in order to obtain a good fit to the data. A similar choice of scales was needed to describe 
double-tagged events at LEP in an analogous NL-BFKL framework [23]. Equations (^) and (|l~6| ) are 



solved using an approximate numerical technique, described in detail in 22] and |1 



4 The exact LL BFKL solution 

The BFKL kernel (|l^) in the leading logarithmic approximation exhibits, in the impact parameter 
representation, invariance under conformal transformations Q. The conformal symmetry of the kernel 
permits the following expansion of the amplitude in the basis of eigenfunctions E n ^ 0]: 

t 2 n ^ f°° , u 2 + n 2 /A 



(2vr) 3 ^ 7_oo [v 2 + (n - l) 2 /4] W 2 + (n + l) 2 /4] 

x exp^^z] /£„(</) (/£,(?))* (19) 

where 

XnH = 47ie(^(l)-^(l/2 + |n|/2 + ii/)J (20) 
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is proportional to the eigenvalues of the BFKL kernel and 

l2, 



f d k f 

In,v{l) = J (2vr)2 Jj4( - fc ' g ' ) J d2 P ld2 P 2 E n , v (pi, P2) exp(ik ' Pi + ~ k ) ■ Pi) ( 21 ) 



(and analogously for the index B). The eigenfunctions are given by 



E n ,u(Pl,P2) 



Pi - Pi 
pi Pi 



Pi - Pi 
P1P2 



*\ h 



(22) 



where h = 1/2 + n/2 + iv and h = 1/2 — n/2 + iv. Here k and g are transverse two dimensional 
momentum vectors, and p\ and p 2 are position space vectors in the standard complex representation 
(e.g. k = k x + iky). The scalar product in this representation is given by, e.g., k- p\ = k*p±/2 + kp*/2. 
The functions I a = an d 1b = &qq are the impact factors (|8|) . 

The quark impact factor in representation ( |2l| ) was found in ||], generalizing the Mueller- Tang 
subtraction Q to non-zero conformal spin; 



/k| 2 \^ /^*\ n/2 r(l/2 + n/2 - ii/) 
TJ \~q) T{l/2 + n/2 + iu) 



(23) 



for even n and = for odd n. 

The impact factor for the 7 — ► V transition is known for n = j|]. We shall generalize this result 
to arbitrary n, which requires evaluating the following integrals 

l2i 



/d k f 
j—^l lV (k,q) J d 2 p 1 d 2 p 2 E ntV {pi,p 2 ) 

x exp(ik*p x /2 + ikp\/2 + i(q* - k*)p 2 /2 + i(q - k)p* 2 /2). 



(24) 



Changing variables to p\ = R + p/2, p 2 = R — p/2 and integrating over d 2 k we get 
C a. 



^ J d 2 pK Q { n \p\) J d 2 RE n , u (R + p/2,R-p/2) exp{iq*R/2 + iqR*/2). (25) 



47T 

The integral over d 2 R was obtained by Navelet and Peschanski [24]. Inserting their result one obtains 
where 



P^TZo — K,u E n>li {p, p )K {q\\\p\) 



-in 



2vr 2 



(26) 



E n ^{p,p*) 



, , 2if / * \ n/2 



r(l - iv + n/2) r(l - *i/ - n/2) 



and 



J n /2-iv(4*P/4) J-n/2-iu(qP*/^} ~ ("If J-n/wtfp/ 4 ) 4/2+^(^74) 

2 4 ^vr 3 r(|n|/2 - + l/2)r(|n|/2 + zV) 



\n\/2 - iv r(|n|/2 + ii/ + l/2)r(|n|/2 - zV) ' 



(27) 



(28) 
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Thus, we have a result in terms of a double integral over d 2 p = \p\ d\p\ d(f). Further, we represent the 
Bessel functions by their power series expansions J a (z) = {z/2) a X]fclo( — 1)* { z /2) 2k /\P(k + l)r(<r + 
k + 1)] and obtain 



Jn/2-iAl*p/4) J-nj1-iv{lP*l^) ~ (-!)" J-n/2+iv(q*P/4) J n /2+iv(qp* /*) 

(~l) k+l (M|p|/8) 2fc+2/ - 2 ^ exp[#(2Z - 2k + n)] 



EE 



T(k + 1) T(l + 1) T(l + 1 + n/2 - iv) T(l + k - n/2 - iv) 



c.c. 



(29) 



L A;=0 Z=0 

Both the 4> and \p\ integrations in eq. (^) are performed term by term in the sums (p9|). All the 
4> dependence of the integrand in ( p5[ ) is due to expression (^) where in the subsequent terms only 
integer powers of exp(i^) appear. Thus, after the angular integration, only terms with n + 2k — 21 = 



contribute, and one of the summations in eq. (2£) may be trivially performed giving 

-2tt 

Jn/2-iu(q*P/^> J-n/2-iv(qP*/4) - (-1)" J- n / 2 +iv(q* p/4) J n /2+iu(qP* /4) 
(-l)N/2 (| g || p |/ 8 )4«+|n|-2^ 



2 'Em 



z=o 



(1 + Z) T(l + l + \n\/2) T(l + I - iv) T(l + I + |n|/2 - ii/) 



c.c. 



(30) 



Note that the odd n contributions are all zero. Using ( |3Q[ ) and performing the remaining |p| integration 
in eq. (|25|), we have 

/>2tt 

|2 



dlpllpl^o^ilpl) 

X J n /2-iM*P/^)J-n/2-iu{qP*/^) ~ (-l) n J-n/2+iv(q*p/4)J n /2+iu{qP*/4) 



4?r ~ (_l)|n|/2 r 2 (3/2 _ ^ + 2 / + | n |/ 2 ) (|g|/ 4g|| ) 4 '+l"|- 2 ^ 

g| T(l + Z) T(l + l + \n\/2) T(l + I - iv) r(l + I + |n|/2 - ii/) 



c.c. 



(31) 



The obtained series is convergent for |q , |/4g|| < 1, where the infinite sum gives the value of the integral. 
We need to continue analytically the result to |(?|/4(/|| > 1. Thus, we represent this infinite sum of 
terms A(l) as an integral over complex I using a Sommerfeld- Watson type transform. A contour C± in 
the complex I plane is introduced enclosing the half-plane lZe(l) > —1/2 where A(l) has no poles. It 
is possible to construct a function T>(1) with a pole structure and with residues such that the contour 
integral of V(l)A(l)/(2iti) along C±, evaluated using the Cauchy theorem, reproduces the initial sum 
over the index Z. It is easy to verify that 

sinvrii/ r(Z)T(l - l)T(l - iv)T(l - I + iv) 



V{1) 



-7T- 



sin7r/ sin7r(Z — iv) 
has the desired properties. Thus one has 



Y(iv)T{l-iv) 



(32) 



47T ~ (_i)H/2 r 2 (3/2 - il/ + 2 l + \ n \/2) (| g |/4g|| )^+l"l~ 2 ^ 



II 1=0 

4ir 



c.c. 



+ Z) r(l + l + \n\/2) T(l + 1 - iv) r(l + I + \n\/2 - iv 

dlV(J) (-I)'"'/ 2 (|g|/4g||) 4/ +l"l- 2w r 2 (3/2-^ + 2Z + |n|/2) 
4 J c 2ni r(l + r(l + Z + \n\/2) T(l + Z - iv) T(l + I + |n|/2 - iv) ' 



(33) 
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The contribution from the region of complex I — > oo vanishes in the limit, so the value of the contour 
integral is given by a line integral 



■ 4?r ~ (_i)H/2 r 2 (3/2 -iv + 2l + \n\/2) (M/^,) 4 ^"!- 2 *" 

- 4 r(1 + r(i + / + M/2) r(i + 1 - w) r(i + 1 + \ n \/2 - iv) 

Air /•— 1/2+ioo j/ 

--(-l) |n|/2 / fr (k|/4g|,) 4i 



C.C. 



9|| J-l/2-ioo 

F(l - iv) T(l - / + iv) Til) T(l - I) T 2 (?>/2 + 2l + \n\/2-iv) 

x r(ii/)r(i - iv) r(i + T) r(i + 1 - iv) r(i + z + \n\/2) r(i + z + h/2 - ^) ' 

This form is suggestive of the n = result of Substituting s = 2Z + 1 — iv, r = q 2 /^q 2 and using 
the Euler T-function relations such as the doubling formula 2 2z ~ 1 T(z)T(z + 1/2) = y/irT(2z) allows 
simplification of the integrand. Inserting ( |27| ) and into (|26|), taking into account (|34|) and using 
the identity (— l)l n l/ 2 r (^+H/2) r(i-^Hn|/2) _ ^ £ Qr even n ^ we obtain the final answer, 

&tt 2 (\q\ 2 X v fq*\ n/2 fl\ lnl/2 T{l/2-iv+\n\/2) 



I OC 



r(l/2 + ii/+ |n|/2) 
ds _i/2+ s +lnl/2 r(l - s - ii/) r(l - a + ii/) 



2?ri r(l - s/2-iv/2)T(l -s/2 + iv/2) 

r 2 (l/2 + s + |n|/2) 

X r(l/2 + s/2 - iv/2 + |n|/2) T(l/2 + s/2 + iv/2 + \n\/2) 

(35) 

for even n and Inf v = for odd n. This agrees with the corresponding expression of Bartels et al. for 
n = Q. The r.h.s. of equation p5| ) is the desired analytic continuation of the sum of the infinite 
power series (c.f. (|3~l|)) which holds for all values of r. Introducing the notation m = n/2, and using 
the result for In,v, we arrive at the amplitude 

i \ H r 7/2 t ™2 



7«=ou /■ ^ v |m| /. 

_F(z,r)=4Ca 2 £ (--) / ^ 



m=— OO 

moo 



(j/2 + ( m _ l/2)2)(z/2 + ( m + 1/2)2) 



ds r i/2+ s +|m| T{1- s-iv)T(l- s + iv 



2-Ki r(l - s/2 - iv/2)T{l - s/2 + iv/2) 

T 2 {l/2 + s + \m\) 



T(l/2 + s/2 - iv/2 + \m\) T(l/2 + s/2 + iv/2 + |m| 

which is equal to 



(36) 



on 2 m =°° f 

?(z,t) = ^ £ ( _i r y 



o m=oo „ 9 9 

(v 2 + (m-l/2) 2 )(v 2 + (m+l/2) 2 ) 

io ° cZs /T\ 1/2+s+lml T(l/2- s/2 + iv/2)r 2 (l/2 + s + \m\)T{l/2- s/2 -iv/2) 
_ ioo 2^i\j) r(l/2 + s/2 - iz^/2 + jmj) T(l/2 + s/2 + iv/2 + \m\) 

Note that the m contribution is equal to the — m contribution. 



(37) 
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Figure 2: The amplitude squared for the diffractive heavy vector meson production off a quark at zero 
rapidity (z = 0): the exact two-gluon exchange result (dots) compared with sums of contributions up 
to conformal spin n max = 2m max (lines) 

5 Properties of the solutions 

The form of the solution to the BFKL equation given by eq. ( |37| ) is simple enough to perform extended 
studies of the amplitude, including the impact of higher conformal spins. The remaining complex 
integrations over s and v are performed numerically. In the analysis we take the real photon case, 
Q2 = o. First, we set z = in which case the LO BFKL amplitude is described by a simple two-gluon 
exchange, 

Hz = ^r)=Cal{^) (38) 

A valuable cross check of our calculation is to investigate how the two-gluon exchange amplitude 
builds up when subsequent higher conformal spin components are being added. In fig. || we show 
curves corresponding to \J-(z,t)/(C a s )\ 2 approximated by partial sums over m in ( |37| ) up to m = 
m max = 0, 1,4 and compare them with the amplitude given by eq. (p8|). It is clear that results coming 
from the two approaches agree. The sum over conformal spins converges quickly to the exact result, 
however more terms are needed for increasing r. Note, that the m = n = component and the exact 
result have about the same absolute values for r 3> 1 but the signs of the amplitudes are opposite. 

Having checked this, we are in a position to study the importance of higher conformal spins at 
z > 0. Thus, in fig. [3| the ratio of J-(z, r) (with m max = 7) to the m = component is plotted for 
various z as a function of r and for various r as a function of z. The relative importance of the higher 
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Figure 3: Comparison of the BFKL amplitude approximated by a partial sum in eq. up to con- 
formal spin \n\ = \2m\ = 14 with the leading conformal spin result, n = 0. Dependencies of the ratio 
are given as a function of a) z and b) r respectively. 
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Figure 4: Contour plot showing the relative error £ of the cross-section obtained in the n = ap- 
proximation to formula $F\) in the (z,t) plane. In this plot £ < 0. The contour labelled da/dt = 
corresponds to values of z and r for which the exact cross-section vanishes. 

conformal spins for the cross-section may be read out from fig. ^. The contours show constant values 
of 

c _ da/dt\ exact - da/dt\ n=0 

da/dt\ n=0 ' W 

giving therefore the relative error of the conventional leading conformal spin approximation in the 
(z, t) plane. As expected, the correction due to m ^ components decreases with increasing z (or 
y) and increases with increasing r (or \t\). Note also the line da/dt = on which the complete 
amplitude changes sign, leading to a dip in the cross-section da/dt\ exSLC t- Such a dip does not appear 
in the leading conformal spin approximation. 

The analytical results may also be used to test the method and approximations used in the nu- 
merical approach to the non- forward BFKL equation ||, [R]]. In fig. |, a comparison between results 
obtained in these two frameworks is given. With good accuracy the numerical solution coincides with 
the analytical one. At larger values of z and r a discrepancy between the two solutions appears. This 
is caused by an upper cut-off Kp = 10 5 GeV imposed on the gluon transverse momenta, which is nec- 
essary for the numerical approach. The very large gluon momenta that are affected by this procedure 
contribute significantly to the amplitudes when the momentum transfer and/or the rapidity is very 
large. 

The non-leading corrections ||] to the BFKL equation have a large impact on the rapidity depen- 
dence of the cross-sections, as already discussed. At high rapidities, this leads to dramatic effects in the 
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magnitude of the cross-section [22]. Therefore it is important to incorporate those non- leading effects 
in the analysis. In this case, the conformal symmetry of the kernel is broken and no exact analytical 
approach is known yet. Fortunately, using the more straightforward numerical method one may ob- 
tain a solution to the BFKL equation beyond the leading logarithmic approximation (BFKL LL+NL), 



given by eq. (|l6|). In fig. g the rapidity y = ln(l/x) (c.f. (12), (16)) and \t\ dependencies of the cross- 
section da/dt (y, \ t\) for J/ip production off a quark are shown. The cross-sections obtained from the 
LL BFKL and BFKL LL+NL are compared. At the leading logarithmic accuracy the value of the 
fixed coupling a s is not constrained and we choose a s = 0.17. With this choice, the LL and LL+NL 
results in the studied window of y and \t\ have a similar overall normalization. In the BFKL LL+NL 
equation (|l6|), we take sq = 0.5 GeV 2 , and the scales of the running coupling according to (0). It 
may be seen in fig. ^a, that the increase for large y is less steep (i.e. the intercept is smaller) for the 
LL+NL case, in spite of taking a rather small value of a s in the LL BFKL equation. To be precise, the 
pomeron intercept is ap = 1.53 for the LL BFKL curves and ap ~ 1.3 for the non-leading solution. 
It is clear, that the non-leading prediction is much closer to the experimental estimates of ap ~ 1.2 
in the vector meson production process. 

Thus, when both the normalization (which also depends on a s ) and the rapidity dependence are 
taken into account, the need for NL effects in the BFKL kernel should become visible. Comparison of 
the shapes in \t\, given in fig. ^3, demonstrates that they are similar in both cases, especially in the 
low \t\ range. The LL+NL curves are steeper because of the running of a s . To conclude, the main 
impact of the non- leading corrections seems to be a reduction of the pomeron intercept. 



6 Comparison with data 

The results of the model calculation described in the previous sections may be compared to the ZEUS 
data H on diffractive J/ip photoproduction. In this measurement, the photon virtuality Q 2 ~ and 
the photon-proton collision energy is in the range 80 GeV < W < 120 GeV. In fig. [7| the data are 
shown together with the theoretical curves obtained from eq. ([!]) with various models. These curves 



were obtained using the CTEQ5L parton distributions [25]. The continuous curve is given by the 
analytical solution of the BFKL equation (|37|) with a s = 0.21 and all conformal spins included, and 
the dashed one corresponds to the leading conformal spin (n = 2m = 0) in (|37|). The dotted curves 
are obtained from the BFKL equation with non- leading effects (|l6[). For the lower curve, the strong 
coupling constant in the impact factors is evaluated at the scales given by (|i~7D, and for the upper 
curve we use the values given by (|T8l) as described in Sec. ||. 

The LL BFKL results fit the data very well and the difference between the leading conformal spin 
and the full solution is small, although the discrepancy increases with t. The non-leading BFKL results 
are in rather good agreement with the data when low scales of a s are chosen (O) , but underestimate 
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Figure 6: The cross-sections do~/dt(y, \t\) for diffractive J/ip production off a quark shown as a function 
of a) rapidity y for \t\ = 3 Ge V 2 (uppermost curves), 10 GeV 2 (in the middle) and 30 GeV 2 (lowermost 
curves), and b) momentum transfer \t\ for y = 4 (lowermost curves), y = 7 (in the middle) andy = 10 
(uppermost curves). Solid curves are obtained from the LL BFKL equation and dashed ones from the 
BFKL LL+NL equation 
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Figure 7: The cross-section da/dt for diffractive J/ip photo-production off proton shown as a function 
of the momentum transfer \t\. ZEUS data are compared with the theoretical results from: LL BFKL 
(continuous line), the leading conformal spin approximation to LL BFKL amplitude (dashed line) and 
the BFKL equation with non-leading corrections (dotted lines where lower and upper lines correspond 
to the choices of scales given by §T\ ) and respectively). A correlated 10% uncertainty of the 
normalization of data points is not included into the error bars. 

the data when the most natural choice ( |17| ) of scales is made. Thus, in formulating predictions 
extrapolating beyond the currently measured kinematical window, we will use the data-guided option 
(p~8|) . We have also checked, that the sensitivity of the normalization of the non- leading BFKL results 
to the value of s$ is relatively large within the considered window: for the scale choice (]17[), it increases 
by a factor of two when sq decreases from 0.5 GeV 2 to 0.1 GeV 2 . The i-dependence is only slightly 
steeper at low t for the latter case. Therefore, the shape of the cross-section is stable against variation 
of details, while the normalization is more uncertain, and critical tests of BFKL should be based on 
the shape rather than the normalization of the cross-section. 

Recall that an important feature of the non-leading BFKL solution is the emerging value of the 
pomeron intercept of about 1.3, to be compared with the LL BFKL value ap = 1.56 for a s = 0.21. 
Thus, we expect that the LL BFKL should overestimate significantly the cross-sections for larger 
average collision energies W. An interesting experimental verification of the impact of non-leading 
corrections on the scattering amplitudes could be provided by performing analogous measurements at 
higher energies W. In fig. || theoretical estimates from LL BFKL and non-leading BFKL are shown 
for the photon-proton collision energy W = 100 and 200 GeV, for a wide t range. The parameters are 
adjusted to give the best fits of the presently available ZEUS data, that is a s = 0.21 for LL BFKL, 
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so = 0.5 GeV 2 and the running coupling being taken at scales (|H|). One may see that although the 
inclusion of non-leading corrections lowers the expected value of the cross-sections at W = 200 GeV 
it may be insufficient to discriminate between the models. We have checked that the impact of higher 
conformal spins may be safely neglected at W = 200 GeV and \t\ < 10 GeV 2 . 

Note that the cross-sections fall off much steeper with increasing t when the non-leading corrections 
are included (see fig. || a,b). This is mostly because of the running of the coupling with the energy 
scale related to the momentum transfer. 

7 Conclusions 

In this paper an analysis was performed of BFKL amplitudes for diffractive heavy vector meson 
photoproduction at large momentum transfer. We obtained an explicit complete solution to the 
leading-logarithmic BFKL equation describing this process. The novel feature of our approach is 
the inclusion of terms subleading at very high rapidity, corresponding to higher conformal spins in 
Lipatov's expansion of the BFKL amplitude. These subleading effects were found to reduce theoretical 
expectations for the cross-sections by about 10% in the kinematical window currently probed by 
experiments on the J ftp production. This result gives a firmer ground for the previous results, based 
on the leading conformal spin approximation. The relative importance of higher conformal spins 
increases, however, with decreasing collision energy or increasing momentum transfer, as shown in 

%• | 

Also non-leading corrections to BFKL equations were taken into account phenomenologically by 
using the running coupling constant and applying the so-called consistency constraint in the BFKL 
kernel. In this case a numerical method was used to solve the equation. The main influence of non- 
leading corrections was found to be a reduction of the hard pomeron intercept to about 1.3, close to 
the value determined from experiment. 

Results obtained from both approaches were compared to the experimental data on the t-dependent 
differential cross-section for J/ifi photoproduction at "fp collision energy W ~ 100 GeV. In both cases a 
good fit was obtained, although we found that the cross-section grows much slower with rapidity when 
non-leading corrections are included, leading to a discrepancy from the LL BFKL results increasing 
with rapidity. The ratio of the differential cross-sections at W = 200 GeV and W = 100 GeV may be 
used to find the influence of non-leading corrections to the BFKL equation if the data are accurate 
enough. 

To summarize, we provide more insight into the BFKL mechanism of diffractive heavy vector 
meson production and confirm that the available data are consistent with BFKL expectations. 
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Figure 8: The cross-sections da/dt for diffractive J/tp photo-production off proton shown as a function 
of the momentum transfer \t\ for ^p collision energy a) W = 100 GeV and b) W = 200 GeV. The 
continuous and dotted lines represent the LL BFKL and BFKL with non-leading corrections results 
respectively. 
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